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INVERSE SPECTRAL PROBLEMS FOR STURM-LIOUVILLE 
OPERATORS WITH SINGULAR POTENTIALS^ 

R. O. HRYNIV AND YA. V. MYKYTYUK 


Abstract. The inverse spectral problem is solved for the class of Sturm-Liouville op¬ 
erators with singular real-valued potentials from the space !)■ The potential 

is recovered via the eigenvalues and the corresponding norming constants. The recon¬ 
struction algorithm is presented and its stability proved. Also, the set of all possible 
spectral data is explicitly described and the isospectral sets are characterized. 


1. Introduction 


The main aim of the present paper is to solve the inverse spectral problem for 
the class of Sturm-Liouville operators with singular real-valued potentials from the 
space IU2"^(0,1). Given a real-valued distribution q G 11^2“^(0,1), we dehne a Sturm- 
Liouville operator T acting in the Hilbert space Ti. := L2(0,l) and corresponding to 
the differential expression 


( 1 . 1 ) 


I : = 




dx"^ 


+ q 


and, say, the Dirichlet boundary conditions by means of the regularization method due 
to Savchuk and Shkalikov [^|. Namely, we take a real-valued cr G 7Y such that 
a' = q in the sense of distributions (thus a is a distributional primitive of q) and put 

(1.2) Tu = T„u = la{u) := —{u' — au)' — au' 


on the domain 


(1.3) T){T„) = {u G lUj^(0,1) \ u' — au e 1), G(m) G H, m(0) = «(!) = 0}. 

Observe that, in the sense of distributions, G(w) = —u" + qu for all u G ©(T^.). In par¬ 
ticular, the operator T„ does not depend on the particular choice of the primitive a and 
for regular (i.e., locally summable) potentials it coincides with the standard Dirichlet 
Sturm-Liouville operator corresponding to O)- Also Tfj depends continuously in the 
uniform resolvent sense on the primitive a E Ti. and thus it is a natural Dirichlet 
Sturm-Liouville operator related with (O) for an arbitrary q = a' E 1U2~^(0,1). Note 
that the class of potentials considered includes, e.g., the Dirac 5-hke and Coulomb 
l/a:-hke potentials that have been extensively used in quantum mechanics and mathe¬ 
matical physics n, 

It is known that for any real-valued a E Ti the operator dehned above is a 
selfadjoint operator with discrete simple spectrum (A^), k E N, and that Xk have the 
asymptotics Xk = Tik + Hk with an f' 2 -sequence (pfc) |^, ^ We recall that for 
regular potentials q the above asymptotics refines to having jXk = 0{l/k). 
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The present paper has arisen as an attempt to answer the following question: Which 
sequences (A^) consisting of real pairwise distinct numbers and obeying the above- 
mentioned asymptotics are indeed spectra of Sturm-Liouville operators with singular 
potentials from hh2~^(0,1)? This has led us to the inverse spectral problem for the op¬ 
erators considered, i.e., to reconstruction of the potential q based on the corresponding 
spectral data. 

In the regular situation knowing only the spectrum (A^) is insufficient: there are 
many different potentials q (called isospectral) producing Sturm-Liouville operators 
with the same Dirichlet spectrum. It was shown by Poschel and Trubowitz [|4 


that the set of all potentials in H with a given admissible spectrum (A^) (i.e., real, 
simple, and obeying the asymptotics \k = Tik + 0(1/k)), is analytically diffeomorphic 
to the weighted space ^ 2 ('R'n) with the weights Wn = n. (A rather exceptional situation, 
where the spectrum determines the potential uniquely, was pointed out by Ambar¬ 
tsumyan [P; namely, he proved that for the Neumann Sturm-Liouville operator on 
(0,1) the equalities Xk = 7r(fc — 1) for all fc G M imply q = 0.) 

Therefore to recover the potential q uniquely some additional information besides 
the spectrum must be supplied. This can be, e.g., knowledge of the potential on half of 
the interval (0,1) 00.001. or the spectrum of a Sturm-Liouville operator given 
by the same differential expression but different boundary conditions [0,J^ , or three 
spectra—one for the whole interval and the others for two halves of it ||11| , |23|| . Another 
kind of information is the squared norms of properly normalized eigenfunctions (the 
so-called norming constants) ||I^ , and that is exactly the setting we shall treat 
in this paper. The aforementioned settings of the inverse spectral problem will be 
considered elsewhere. 

A complete solution of the inverse spectral problem for a class of Sturm-Liouville 
operators must consist of two parts: (1) an explicit description of the set SV of spectral 
data for the operators in SC and (2) development and justihcation of the method of 
recovering the operator in SC that corresponds to arbitrary given spectral data in SV. 

The algorithm of recovering the potential q from the spectral data of a regular 
Sturm-Liouville operator based on the transformation operators and the so-called 
Gelfand-Levitan-Marchenko equation was developed by Gelfand and Levitan Q 
and Marchenko in early 1950-ies. In particular, a reconstruction method via 
the spectrum (A^) and the norming constants (afc) was suggested in ^ (see also |1l9|), 
though no precise description was given therein of the set of all spectral data generated 
by the considered class of regular Sturm-Liouville operators. An alternative method 
for reconstruction of the potential q by two spectra was developed by Krein |]^. A 
different approach was suggested later by Poschel and Trubowitz for the class of 
Sturm-Liouville operators with potentials from Ti. in [Q . The authors studied in detail 
the mapping between the potentials in H and the spectral data, proved solvability of 
the inverse spectral problem and, in particular, completely characterized the spectral 
data. 

In 1967, Zhikov considered the singular case where the potential q is the de¬ 
rivative of a function of bounded variation (i.e., where g is a signed measure). The 
corresponding Sturm-Liouville operator was dehned through the equivalent integral 
equation, and complete solution to the inverse spectral problem (in particular, neces¬ 
sary and sufficient conditions on spectral data) was given. In fact, in the problem 
was reduced to recovering q by the corresponding spectral function. 

Singularities of different kind were treated by Rundell and Sacks in [^. They 
considered a Sturm-Liouville operator in impedance form Su = -{au')' in the space 
L2((0, l);a) with a positive impedance a G W^lO, 1]. With the help of transformation 























operators they found necessary conditions on the spectral data, solved the inverse 
spectral problem, and suggested a numerical algorithm for finding a solution. Observe 
that S is similar to Tu = —u" + qu with q = {y/a)”/-/a. In particular, for a G 1] 

we get q G though there are many such functions a producing the same 

qeW,-\o,i). 

Coleman and McLaughlin |P] took a different approach to the inverse spectral 
problem for Sturm-Liouville operators in impedance form. Namely, they recast the 
equation {au')' + Xau = 0 in the form u" + bu' + \u = Q with b := a'/a and then 
modihed the method of accordingly. In particular, they described the set SV of 
the spectral data and proved unique solvability of the inverse spectral problem for the 
case where b G ^ 2 ( 0 ,1). 

Note that other types of singularities (e.g. discontinuous a for impedance Sturm- 
Liouville operators S, 1/x'^-like potentials) were treated by Andersson 0, Carl- 
Hald 


SON 


14 , McLaughlin 


Yurko 29 


a.o. 


Here, we generalize the classical approach due to Gelfand, Levitan, and Mar¬ 
chenko and completely solve the inverse spectral problem for the class SC of Sturm- 
Liouville operators with singular potentials from 1 ^ 2 "^ (0,1). Namely, we give an explicit 
description of the set SV of spectral data, explain how to recover q from an arbitrary 
element of SV, and study dependence of q on spectral data. As a byproduct, we show 
that the set of potentials in 1 ^ 2 "^ (0,1) that are isospectral to a given one is analytically 
diffeomorphic to the space £ 2 - The main tool of the reconstruction procedure is the 
transformation operators for the class of singular Sturm-Liouville operators in SC that 
were constructed in . 

Our primary purpose in this paper is to treat the general singular case. We do not 
prove here that if the spectral data have better asymptotics, then the recovered poten¬ 
tial is smoother—e.g., if the spectral data formally correspond to a regular potential 
from Ti, then the recovered potential indeed falls into Ti. In a subsequent work we shall 
justify a more general claim that the reconstruction algorithm suggested here solves 
the inverse spectral problem for the class of Sturm-Liouville operators with potentials 
from 11 ^ 2 "“^(0,1) for every fixed a G [0,1]. 

The organization of the paper is as follows. In the next section we exploit transfor¬ 
mation operators to study the direct spectral problem, i.e., to give necessary conditions 
on the set SV of spectral data for the class SC of Sturm-Liouville operators considered. 
In Section ^ we establish connection between the spectral data and the transformation 
operators and, in particular, derive the so-called Gelfand-Levitan-Marchenko (GLM) 
equation. In Section || the GLM equation is proved to possess a unique solution for 
any element from SV, and this solution is shown to be the transformation operator for 
some Sturm-Liouville operator from 5/1 in Section Moreover, the element from SV 
that we have started with turns out to be the spectral data for the Sturm-Liouville op¬ 
erator found, so that the inverse spectral problem is completely solved. In Section ^ we 
show stability of the reconstruction algorithm and characterize isospectral sets, and in 
the last section we comment on changes to be made for the case of Dirichlet-Neumann 
or third type boundary conditions. Finally, two appendices contain necessary facts on 
Riesz bases of sines and cosines in functional Hilbert spaces and on Hilbert-Schmidt 
operators respectively. 

Throughout the paper, a G 7/ will denote a real-valued distributional primitive of the 
potential q G 1 ^ 2 ’“^(0,1), will stand for the quasi-derivative u' — au of a function u, 
and Hull will denote the 7/-norm of u. 
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2. Transformation operators and direct spectral problem 


In this section we shall solve the direct spectral problem for Dirichlet Stnrm-Lionville 
operators T^. with a G i. e., we shall describe the set SV of spectral data for T^.—the 
seqnences of eigenvalues (A^) and norming constants (a^) introduced below—when a 
runs over Ti. The main tool will be the transformation operators. 

Suppose that a G is real-valued; we denote by a Sturm-Liouville operator 
— ^ with the Dirichlet boundary condition at the point a; = 0. More precisely, 
To- acts according to 

To-m = la{u) = — — ctmW — a'^u 

on the domain 


D(To) = G W^[0, 1] I mW G W^[0, 1], IM G T2(0, 1),m(0) = 0}. 

(We recall that := u' — au is the quasi-derivative of u.) According to the dehnition 
of /o, the equation la{u) = n is understood in the sense that 


d (u\ / a 

dt VmW ) \^-cr^ 





1 


so that its solutions enjoy the standard uniqueness properties. 

One of the main results established in [|^ is that the operators T„ and Tq possess 
the transformation operators, which perform their similarity. 


Theorem 2.1 (0]). Suppose that a G 7i; 
operator K„ of the form 


then there exists an integral Hilbert-Schmidt 


( 2 . 1 ) 


(K„u)(x)= k,(x,y)n(y) dy 

Jo 


such that I is the transformation operator for To- and Tq, i.e., 


( 2 . 2 ) 


TAI + K,) = (/ + 


If a is real-valued, then the kernel fco- of is real-valued too. For every x G [0,1] and 
y G [0,1] the functions ka{x, •) and ka{-,y) belong to H and the mappings x h-> k„{x, •) 
and y i—> kcj{-,y) are continuous in the U-norm. Moreover, the operator with 
properties (0-(iD is unigue. 


Remark 2.2. Some other properties of the transformation operator established 
in imply that I -|- K„ preserves the initial conditions at x = 0 in the sense that for 
any u G 1^2^[0,1] and v ■.= {I -\- K„)u we have n(0) = m(0) and nW(0) = m'(0). 


i.e. 


Denote by T^ the restriction of T„ by the Dirichlet boundary condition at a; = 1, 
To- := To-|j)o, where Dq = {u & Ti(To-) | m( 1) = 0}. It is known |2^ 


that To- is 

a bounded below selfadjoint operator with simple discrete spectrum. Denote by A^, 
A; G N, eigenvalues of To- in increasing order. We may assume that all the eigenvalues 
are positive (and that such are A^), as otherwise this situation can be achieved by 
adding a suitable constant to q. 

Let also Uk, k E N, he the eigenfunction of To- corresponding to the eigenvalue A^ 
and normalized by the condition ^^^(0) = V^Xk] then Ok := is the corresponding 

norming constant. Notice that in view of Remark we have Uk = {I -h Kcr)vk, where 
Vk{x) = A/2sinAfcX. We also put Vkp{x) = \/2sm.TTkx. 

The following statement was proved in |^, | 
the sake of completeness. 


281, but we sketch its proof here for 
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Lemma 2.3. Suppose that a Ehi is real-valued and that < ... are eigenvalues 

ofT^ with \1 > 0. Then A^ = vr/c + fik, where the sequence {pik)'^=i belongs to £ 2 - 


Proof. Since Uk{x) = Vk{x) + k^{x, y)xk{y) dy, the numbers Xk are zeros of the func¬ 

tion 


<h(A) := sin A-I- / y) sin Xydy. 


Recall that k^{l, ■) E TC; whence fccr(l, y) sin Xydy ^ 0 as A —*■ cx) by the Riemann 
lemma. Rouche’s theorem then gives Xk = 7rk-h p,k with /i^ —0 as A; —> cxo (see details 
in 1^ Ch. 1.3] and |0). Therefore (see Appendix a {sin Afcx}^^ is a Riesz basis 
of Pi and the numbers 


sin Afc = (-l)^sinpfc = 


ka{l,y) sin Xkydy 


are the Fourier coefficients of — A:o-(l, ^ ^ fhe biorthogonal basis. It follows that 
(sinpfc) G £ 2 , hence (/ifc) G £2 and the proof is complete. □ 


The next lemma gives the asymptotics of the norming constants 


Lemma 2.4. Suppose that a G PC; then the norming constants ak are of the form 
1 + j3k, where the sequence {(3k) belongs to £ 2 - 

Proof. Observe that the vectors Ufc, k G N, form a Riesz basis of PL since so do Vk 
and Uk = {I + Kfj)vk with bounded and boundedly invertible / -|- K^. Therefore the 
Ti-norms of Uk are uniformly bounded, and in view of the relations 

\(3k\ = Ill^fclP - |kfc,o|P| < (1 + l|wfc||)l|wfc - 

the lemma will be proved as soon as we show that 

CXD 

- VkpW^ < 00 . 

k=l 

This inequality follows from the representation 

Uk - Vk,o = {k + K^){vk - Vk,o) K^Vkft 

and the fact that both sequences (||(/-|- K„){vk — 'yfc,o)||) and (||iFo-nA:,o||) belong to £ 2 - 
the former due to the relation 

Vk{x) -Vkp{x) = 2A/2sin YCOs(7rfc-h/ifc/2) = 0(|pfc|) 

and the inclusion (/ifc) G £2 (see the previous lemma), and the latter because is a 
Hilbert-Schmidt operator and (n^.o) is an orthonormal basis of PL (see Appendix |^. 
The proof is complete. □ 

Denote by SV the set of all pairs of sequences {(A^)^]^, (a^)^^} satisfying the 
following two conditions: 

(Al) Xk are all positive, strictly increase with k, and obey the asymptotic relation 
Xk = Tik + pfc with some ^ 2 -sequence {yk)T=P 
(A2) ak = I + (3k > 0 with some ^ 2 -sequence (/3fc)^i. 

Also SC will stand for the set of all positive Dirichlet Sturm-Liouville operators T^j 
with a E PL real-valued. Lemmata |2.3| and |2.4| demonstrate that the spectral data for 
the operators in SC belong to ST). Our next task is to show that, conversely, any 
element of ST) is spectral data for some operator in SC. 
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3. Connection between spectral data and transformation operators 


In this section we shall derive a relation between the spectral data for a Dirichlet 
Sturm-Liouville operator T„ with a & Ti and the transformation operator K^. This 
relation will be used in the next section to find given the spectral data and thus to 
solve the inverse spectral problem. 

Suppose therefore that {(A^), (cifc)} G SV is the spectral data for a Sturm-Liouville 
operator with a E H. We recall that Afc = vr/c -|- /i^ > 0 and ak = I + Pk > 0 with 
some ^ 2 -sequences (/ifc) and (Pk)- Put 


N 


U :=s-hmy^— {■,Vk)vk, 

Oik 

k=l 


where Vk{x) = \/2sm.XkX and s-lim stands for the limit in the strong operator topology 
of 7i. Since by Proposition |A.1| the system {vk)^=i forms a Riesz basis of Ti and 
inffcgN l/ofc > 0, the operator U is bounded, selfadjoint, and uniformly positive. 


Lemma 3.1. For all j,k G N, the following relation holds: {U ^Vj,Vk) = (Xkdjk, where 
6jk is the Kronecker delta. 


Proof. Note that 

N 

U~^ = s-lim akp, Wk)wk, 

k=l 

where (wk) is a basis biorthogonal to {vk) (see 
Therefore 

N 

{U~^Vj,Vk) = s-\imy^ai{vj,wi){wi,Vk) = aidjidik = akdjk, 

N^OO ^ 

1 = 1 

and the proof is complete. □ 


12, Ch. VI] and also Appendix 


Lemma 3.2. F := U — I is an integral operator of the Hilbert-Schmidt class with 
kernel 

(3.1) f{x,y) = (j){x+ y) - (j){x-y), 
where 

(3.2) 0(s) = ^^^costt/cs -cosAfcsj 

is an L 2 { 0 , 2)-function. 

Proof. Recall that we have denoted by Vkp, k E N, the function a/2 sin tt/ ex. Since 
I = s-limAT^oo Y.k=ii'Hk,o)vk,o, we have 

U - I = s-limy^(— {■,Vk)vk - {■,Vk,o)vkft). 

k=l 

Observe that — {■,Vk)vk — {■,Vko)vko is an integral operator with kernel 
Oik 


fk{x,y) 


— sin XkX sin Xky — 2 sin nkx sin irky 
Oik 

= cos7r/c(x + y) -cos Afc(x + y) — cos7r/c(x — y) -\ -cos Afc(x — y). 

Oik Oik 
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Letting 


N 


07 V := ^^(cosvrfcs-cosAfcsj, 


we see that 


k=l 


N 


Fn ■= y^( - {■,Vk)Vk - {■,Vk^o)Vk,0 

is an integral operator with kernel fN{x,y) := (j)N{x + y) — — y). 

Put Pk := Pk/oik] then 1/ak = I - Pk with (0fc)^i e £2 and 

cosnks -cos AfeS = 2 sin(p,fc'S/2) sin[( 7 rA; + /ifc/2)s] + Pk cos A^s. 

Oik 

Recall (see Corollary |A.2|) that both (sin[( 7 r/c+ /i7./2)s])^^ and (cosA^s)^]^ constitute 
Riesz bases of their closed linear spans in L 2 ( 0 , 2 ). Since, moreover, 

2 sin(/ifcs/ 2 ) = s/ifc + 0{\yk\^) 


as k ^ 00 , the series for 0 is convergent in L 2 ( 0 , 2 ) and hence pN —^ 0 € L 2 ( 0 , 2 ). 
Thus /at converge in the -L 2 (( 0 ,1) x (0, l))-norm to / given by (|3.1|) . This implies that 
Fjsf converge in the Hilbert-Schmidt norm to the integral operator F with kernel /, 
and the lemma is proved. □ 


Lemma 3.3. (/ + K^){I + F){I + K*) = F 

Proof. Since {uk /is an orthonormal basis of Ft and Uk= {I + Ka)vk, we have 


TV 1 TV ^ 

I = s-lim"^ — {■,Uk)uk = s-limy^ —(•, (J + K„)vk){I + K^)vk 

TV—>00 ^ Ctk TV—>00 Oik 

k=l k=l 


N 


= {I + K^) s-hmV—(•,n,K {I + K:) = {I + K,){I + F){I + K:). 

Ltv^oo ak J 

k=l 

The lemma is proved. □ 


Corollary 3.4. The kernels ka{x,y) and f{x,y) of the operators K„ and F satisfy the 
following Gelfand-Levitan-Marchenko (GLM) equation for a.e. x,y with x > y: 

(3.3) ka{x,y) F f{x,y) + j k^{x, s)f{s,y) ds = t). 

Jo 

Proof. Observe that Kf is an integral Volterra operator with upper-diagonal kernel 
K{x,y) = ka{y,x), i.e., 

{K*^u){x)=f kl{x,y)u{y)dy 

J X 

and kf{x, y) = 0 for x > y. Therefore the operator I + K* is invertible and its inverse 
can be written in the form I + K, where K is an integral Volterra operator with 
upper-diagonal kernel. By Lemma ^.3| , 

{I + K,){I + F) = {I + Kl)-^ = I + K- 

spelling out this equality in terms of the kernels k^ and / for x > y, we arrive at the 
GLM equation ( p.3|) . □ 
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4. Solution of the GLM equation 


In the previous section, we showed that the spectral data {(A^), (ak)} for the Sturm- 
Liouville operator are connected with the transformation operator for the pair 
Tfj and Tq through the GLM equation ( p.3|) . In this section we shall prove that the 
GLM equation is uniquely soluble for K^. 

Recall (see Appendix P) that the ideal ©2 of Hilbert-Schmidt operators consists of 
integral operators with square summable kernels and that the scalar product {X, Y )2 '■= 
tr^XY*) introduces a Hilbert space topology on © 2 . Denote by ©^ the subspace of 
©2 consisting of all Hilbert-Schmidt operators with lower-triangular kernels. In other 
words, T G ©2 belongs to ©^ if the kernel t of T satishes t{x, y) = 0 for 0 < x < y < 1. 
For an arbitrary T G ©2 with kernel t the cut-off L*" of t given by 


t+{x,y) 


t{x, y) foT X >y 

0 foT X < y 


generates an operator G ©2 , and the corresponding mapping : T 1 —>■ turns 

out to be an orthoprojector in ©2 onto ©^, i.e. and {V^X,Y )2 = 

{X,V^Y )2 for any X,Y & © 2 ; see details in |]^, Ch. 1.10]. 


With these notations, the GLM equation ( p.3|) can be recast as 
(4.1) K + V^F + V^{KF) = 0 


or 

(J + V^)K = -P+F, 

where Vx is the linear operator in ©2 dehned by VxY = V'^{YX) and I is the identity 
operator in © 2 . Therefore solvability of the GLM equation is closely connected with 
the properties of the operator Vp. 

Lemma 4.1. The operator Vx depends continuously on X E ^{TL). 

Proof. This is a straightforward consequence of the fact that Vx depends linearly on X 

n Gi "t” n "f" 

||F+r||e, <rX||e, <||X||||r||e„ 

see 1^, Ch. 3]. □ 

Lemma 4.2. Suppose that I + X is a uniformly positive operator in TL. Then I -|- Vx 
is uniformly positive in 

Proof. Given any Y G ©^, we observe that 

(F + F+y, F )2 = (F, F )2 + (FX, F )2 = tr (F(/ + X)F*). 

If e > 0 is chosen so that I + X > el, then Y{I + X)F* > sYY* and by monotonicity 
of the trace we get 

{{I + V^)Y,Y)^>e{Y,Y)^. 

The lemma is proved. □ 

Lemma 4.3. Suppose that sequences and {ak)'^=i satisfy conditions (Al) and 

(A2) and that F is an integral operator with kernel f of (0). Then equation (H) 
has a unique solution K, which belongs to ©^. 

Proof. Under assumptions (Al) and (A2) the operator F is of Hilbert-Schmidt class 
and 





(cf. the proof of Lemma Moreover, U := J + F is uniformly positive as noted at 
the beginning of Section ^Therefore the operator X + Vp is uniformly positive in 
by Lemma and 

K := -(X + X+)-^X+F G 6+ 


is a unique solution of (|4.1|). The lemma is proved. 


□ 


Remark 4.4. It is well known that a necessary and sufficient condition for the GLM 
equation to he soluble is that, for any a G [0,1], the operator 

u^u{x)+ f{x,y)u{y)dy 
Jo 

is houndedly invertible in H. For a symmetric kernel f this condition is easily seen to 
be equivalent to uniform positivity of the operator I + F. 

Also, the GLM equation is uniquely soluble if and only if I + F can be factorized as 
(/ + X+)(J + X“) with some X~^ G 6^ and X~ G ©2 see |T^, Ch. IV]. 


5. The inverse spectral problem 


In the previous two sections we showed that the operator F constructed via the spec¬ 
tral data {(A^), (a^)} of a Dirichlet Sturm-Liouville operator T„, a E H, is connected 
with the transformation operator I + K„ for and Tq through the GLM equation (|3.3|) 
or (|4.1|) and then proved that the GLM equation is uniquely soluble for K^. Suppose 


that instead we have started with an arbitrary element {(A^), («*;)} G SV, constructed 
F in the same manner, and found a solution K of the GLM equation; is then I + K a. 
transformation operator for some T„, a E U, and Tq? 

In this section we give an affirmative answer to this question and hnd an explicit 
formula for the corresponding a. Moreover, we show that the element {(A^), (a^)} 
we have started with is indeed the spectral data for the Sturm-Liouville operator Tg. 
with the a found. This completes the solution of the inverse spectral problem for the 
considered class of Sturm-Liouville operators. 

Theorem 5.1. Suppose that {(A^), (oa:)} is an arbitrary element of SV and that F is 
an integral operator with kernel f of O- Let also K be a (unique) solution of the 
GLM equation (^T|). Then there exists a (unique up to an additive constant) function 
a eU such that I + K coincides with the transformation operator I + for the pair 
of Sturm-Liouville operators T„ and Tq. 

Proof. We shall approximate F in the © 2 -norm by a sequence {Fn)')(Li of Hilbert- 
Schmidt operators with smooth (say, inhnitely differentiable) kernels /„ so that the 
following holds: 

(a) the solutions Kn of ( [4.1|) for F replaced with Fn converge to iL in ©2 as n —>• 00 ; 

(b) for each n G N there exists an eU such that I+Kn is a transformation operator 
for the pair T^.^ and Tq; 

(c) {an)((Li is a Gauchy sequence in H. 


Put a := imvan'i then by the operators I + Kn converge in Ti to an operator I-\-K„, 
which is the transformation operator for the pair T„ and Tq. Thus K = K„ yielding 
the result. The uniqueness of a up to an additive constant is obvious. 

The details are as follows. We put 


fn{x, y) := (fnix + y)- (t)n{x - y), 
9 










where 


0n(-S) := ( cos Tlks -cos AfcS ). 

We recall (see the proof of Lemma p.2|) that 0^, —0 in ?-f as n —oo. This implies 
that fn converge to / in -L 2 (( 0 ,1) x (0,1)), i.e., that converge to F in ©2 as n — 
cx). The corresponding operators / + are uniformly positive so that the equations 
Kn + Fn + V^{KnFn) = 0 have unique solutions Kn = —(X + Vp^)~^V^Fn G 
Combining all these statements, we hnd that 

||(I + P+)-‘-(J + Pp-‘||s,s,)-0 


and 


||A'„ - A'lls, = 11(1 + V+J-^V*F„ - (I + P+)-'P+A||e. 

< 11(1+ ?’?.)■' - (I + 'Pj)-‘lls(e.)ll'P*A„||s, 

+ 11(1 + Pj!:)-‘||B(e.)ll'P+A„ - V*F\\e, ^ 0 

as n —> CX) thus establishing (a). 

Note that the kernel of Kn solves the GLM equation 

kn{x,y) + fn{x,y)+ kn{x, s)fn{s,y) ds = 0, x>y. 

Jo 

Since fn is smooth, the classical result []^ Ch. II], [^, Ch. 2.3] states that kn is at least 
once continuously differentiable and that I + Kn is the transformation operator for the 
Sturm-Liouville operators Ta-„ and Tq, where cr„ is a primitive of qn{x) := 2-^kn{x,x). 
The suitable choice of the primitive (suggested by the GLM equation above) is 

( 5 . 1 ) an{x)-.= 2kn{x,x)+ 2(j)n{0) =-2(j)n{2x)-2 kn{x,s)fn{s,x)ds, 

Jo 

and (b) is fulhlled. 

To prove (c), we observe that 

II- CTmir < 12 f Ifni^x) - (()m{2x)\‘^ dx 


+ 12 dx 


+ 12 dx 


{knix,y) - kmix,y))fniy,x) dy 
kmix, y) {fn{y, x) - fm{y, x)) dy 


Since fn{,x, y) = (/)„(x + ?/) — fn^x — y) and (fn form a Gauchy sequence in ^ 2 ( 0 , 2), we 
hnd that 


\f„(y,x)\‘dy <2 \4>„(y + x)\‘^dy + 2 \4>n(.y - x)f dy 


(5.2) 


r-2x 


= 2 


|0n(s)pds<2/ |0n(s)pds<G 
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for some positive C independent of n G N and x G [0,1], and hence 
{K{x,s) - km{x,s))fn{s,x)ds^ 

< / dxH \kn{x,s) - km{x,s)f ds \fn{s,x)\‘^ ds^ 


10 '-JO 
>1 rl 


< c 



0 JO 


\kn{x, s) — kfn{x, s) 1^ ds dx 


= C \\ Kn - K ,^\\%^^0 

as n, m —>■ oo. In the same manner it can be shown that 


Ifni^s^x') /j 7 i(s,x)| ds ^‘2 1 |0^(s) 0m('5)| ds 2||0^ 0 ™|Il 2 (O, 2 ) 


and 


/ k^{x,s){fn{s,x) - f^{s,x)) ds^ < 2||i^mi||j|0n - 0 m||i 2 (O, 2 ) ^ 0 

as n,m —> oo. Combining the above relations, we conclude that {an) is a Cauchy 
sequence in Ti. Therefore (c) is satished and the proof is complete. □ 


Remark 5.2. The arguments used to prove (c) above also justify passage to the limit 
in the TL-sense in (O), and this results in the eguality 

PX 

(5.3) a{x) =—2(j){2x) — 2 / k{x, s)f{s, x) ds, 

Jo 

where k is the kernel of K. If k and 0 are smooth, then the GLM eguation implies that 
a(x) = 2k(x,x) — 20 ( 0 ), thus yielding the classical relation 


for the potential q. 


q{x) = 2—k{x, x) 
dx 


To complete the arguments, we have to show that the spectral data for T„ with 
a eTL just found coincide with the data {(A^), (ofe)} G SV that we have started with. 

Theorem 5.3. {(A^), (ofc)} are the spectral data for T„. 

Proof. Put Wk := (/ + K)vk] then tCfc(O) = 0, ^^'(0) = V2Xk, and, due to similarity of 
To- and Tq, T^jWk = X\wk. The system {wk)'^=i is complete in TL since such is (ufc)^! 
and / + iC is a homeomorphism. Also, Lemmata and imply the orthogonality 
relation 

{wj,Wk) = ((/ + K*){I + K)vj,Vk) = {U~^Vj,Vk) = akSjk, 
in particular, ak = HwfciP- 

It remains to show that Wk are eigenfunctions of T^-, i.e., that Wk{l) = 0. Observe 
that 

(5.4) 0 = {f^Wj,Wk) - {wj,f^Wk) = -M;j^^(l)M;fc(l) + M;j(l)M;fc^(l). 

If Wj{l) = 0 for some j G N, then wj^'(l) ^ 0 due to uniqueness of solutions of the 
equation la{y) = X'jy, and the above relation shows that Wk{l) = 0 for all A; G N as 
required. 

Otherwise Wj{l) ^ 0 for all j G N, and equation ( |5.4|) implies that there exists a 

constant h G M such that r(;j^'(l)/tCj(l) = h for all j G M. Thus A^ are eigenvalues 
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of the Sturm-Liouville operator which is the restriction of by the boundary 
condition = hy{l). However, the eigenvalues of are known to obey the 

asymptotics Xu = 7r(fc — 1/2) + o(l) as k ^ oo (see Section 0), and this contradiction 
eliminates the possibility that tCj(l) do not vanish. 

Thus we have proved that every function Wk satishes the Dirichlet condition at 
X = 1 and so is an eigenfunction of the Sturm-Liouville operator corresponding to 
the eigenvalue A^. Moreover, has no other eigenvalues since the system {wk)^=i 
is already complete in Ti. Henceforth the element {(A^), (afc)} G SV is indeed the 
spectral data for the Sturm-Liouville operator with a of ( |5.3|) , and the theorem is 
proved. □ 

To sum up, we have established the following: 

Corollary 5.4. Two sequences (A^) and (a^) the spectral data for a (positive) 
Sturm-Liouville operator with potential q = a' from the space 1) if and only 

if assumptions (Al) and (A2) are satisfied. 

The corresponding operator T„ is uniquely recovered from the spectral data through 
formula ( b.3|) , in which the functions 0 and f are given by (0 and (O respectively, 
and k is the kernel of the solution K of the GLM equation m- 


6. Stability and isospectral sets 

In this section, we would like to study the correspondence between the Dirichlet 
Sturm-Liouville operators with real-valued a E Td and their spectral data in more 
detail. Namely, we shall show that the potential q = a' E hL2~^(0,l) (and thus the 
operator T^- G SC) depends continuously on the spectral data {(A^), (ofc)} G SV and 
that the isospectral sets are analytically diffeomorphic to the Hilbert space £ 2 - Similar 
results were established in ^ for the regular case q eTL, and in [Q for impedance 
Sturm-Liouville operators. 

First we have to introduce the topology on the set SV. Recall that by dehnition any 
element {(A^), (a^)} G SV is uniquely determined by two ^ 2 -sequences (/r^) and (0^) 
through the relations \k = T^k + yk and ak = 1 -f- 0^. Therefore we can identify SV 
with an open subset of the space £2 x £2 in the standard coordinate system ((/Xfc), (/3fc)) • 
In this way SV becomes a subset of a Hilbert space and inherits the topology of that 
space. 

We shall study the correspondence between spectral data {(A^), ( 0 ;^)} G SV and 
operators To- G SC through the chain 

I—»• 0 I—^ R" I—^ O' e-*• To-, 


in which 0 G T2(0, 2) is the function of (p.2|) and iF G ©2 is the operator of ([4.1|) . 


Lemma 6.1. The function 


0(s) = 

fceN 


(cosnks - 

V ak 



depends continuously in L2{0,2) on the spectral data {(A^), ( 0 ^)} G SV. 
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Proof. We can rewrite the function 0 in terms of the f' 2 -sequences (/ifc) and {(5k) with 
(5k-.= l- l/ak = 13k/ak as follows: 

0(s) = ^^(cosvrfcs — cos(7rfcs + Hks) + (3k cos(7rfcs + /ifcs)) 
fceN 

= 2 sin(/ifcs/2) sin[( 7 rA; + /ifc/2)s] + (3k cos{7iks + /ifc-s). 

fceN fceN 


Therefore it suffices to prove that the mappings from £2 x £2 into -^ 2 ( 0 , 2) given by 


{{^k), (bfc)} ^ V’i(5) := ^7fcCOs(7rfcs + Uks) 

km 

and 

{(t'fc), (7fc)} V’2(s) := ^7fcSin(7rfcs + z/^s) 

km 

are continuous. 

Suppose that {{kk), ( 7 fc)} is another element of £2 x £2 giving rise to a function 
•= SfcGN cos(7rfcs + z/fcs) e 1 . 2 ( 0 , 2); then 


11 V’l -'01 11^2(0,2) < II - %) cos(7rfcs + r'fcs)||L 2 (o, 2 ) 

( 6 . 1 ) 

+ II 5^7fc[cos(7rfcs + z/fcs) - cos{7iks + hs)]\\ ^2(0,2)- 

We may assume that the sequence {irk + Uk) strictly increases; then (cos( 7 r/cs + z/fcs))^^ 
is a Riesz basic sequence in L 2 ( 0 , 2 ), and the hrst summand above is bounded by 
^ll(7fc) — ( 7 fc)l 02 with C = C(^{h'k)) being the corresponding Riesz constant (see Ap¬ 
pendix ^). In view of the inequality 


cos{7iks + Uks) — cos{7iks + z>fcs)| < Ik'k — kk\, s G ( 0 , 2 ), 


the second summand in (|6.1|) is bounded by ||( 7 ^) 102 II(^fc) — {kk)\\e 2 : and the statement 
for 01 is proved. Continuity of the second mapping is proved analogously. □ 


Lemma 6.2. The solution K & &2 of the GLM equation ( [4.1|) depends locally analyt¬ 
ically on <p ^ 1 . 2 ( 0 , 2). 

Proof. The operator F depends linearly on 0 G 7 . 2 ( 0 , 2) and, in view of inequality (p.2| ), 
||-P"|| < 11-^1162 < 2II011 1 ^ 2 ( 0 , 2 ), so that F is an analytic function of 0 (see ||, Ch. 2] 
or 1^ App. A] on analytic mappings of Banach spaces). Next, by Lemma the 
operator Vp G 23 ( 62 ) is continuous in F G 62 (and thus analytic in view of linearity), 
while the inverse function (X — Vp)~^ is locally analytic in Vp. Combining these 
statements and recalling the formula 

R = -(X + X+)-iX+F, 

we easily derive the claim. □ 


Lemma 6.3. The function a of (|5.3|) depends locally analytically inTi on(f) in X2(0, 2). 

Proof. By dehnition a{x) = —20(2a;) — 2 k{x, s)f{s, x) ds, where k is the kernel of K 
and /(s, x) = 0(s + x) — 0(s — x). The second summand above is a continuous bilinear 
function of K and 0 (see the proof of Theorem ^.1|) . Therefore a is a jointly analytic 
function of 0 and K p, Ch. 2], and in view of Lemma ^ the result follows. □ 

13 





We denote by S+ the set of all real-valued a E H, for which the operators T„ are 
positive. Observe that the operators T^., and the spectral data do not change if 
a is replaced by a -|- c with any real c, so that these objects depend in fact on an 
equivalence class a in S^/M (i.e., on the common derivative q = a' for a E a) rather 
than on a G It is proved in |^, ^ that this dependence is continuous for T„ and 
Kfj] we shall show next that also the spectral data depend continuously on d G S^/M. 


Lemma 6.4 (cf. Lemma 4.3]). The mapping 3 d i—>■ (A^) is continuous. 

Proof. We recall that Xk = T^k + pik for some ^ 2 -sequence (pfc) and that the continuity 
of (A^) is understood as continuity of (pfc) in the ^ 2 -topology. Also, A^ are simple zeros 
of the function ^ 

<h(A) = <h(A, a) := sin A-|- / ka{l,y) sin Xy dy, 

Jo 

where fco- is the kernel of the corresponding transformation operator K^r. 

It can be shown (see also m Ch. 1.3]) that 


|<h(A) - sin A| = / k„{l,y) sin Xydy 

Jo 

<I)'(A) - cosA| = / yk^{l, y) cos Xydy 


0 , 


0 


as A —> cx) inside the strip {z E C | |Imz| < 1}. Therefore there exist £o > 0 and a 
sequence (Cn) of circles := {z G C | \z — A„| = r„}, G (0,1/2), such that the 
discs Dn :={. 2 GC| \z — Xn\ < rji} are pairwise disjoint and for all n G N we have 

( 6 . 2 ) 


min |<I)(A)| > £o, 


min |$'(A)| > £o- 
agd„ 


It is proved in that the mapping E+ 3 a h->• fco-(l, •) G TC is continuous, whence 
for any e: G (0, £o) there exists 5 > 0 such that for every d G S’*' from a ^-neighbourhood 
of a we have 11^5(1, •) — k„{l, •)|| < e/4. We £x such e and a function d, denote by 
A^ the corresponding eigenvalues and put ‘h(A) := <I>(A,d) and g := ka{l, •) — ^^-(l, •). 
Simple computation shows that 


(6.3) 

(6.4) 


min |*(A) - *(A)| 

AGOn 



\g{y)\ dy < 60 / 2 , 


min |<hhA) 
AeDn 


<h'(A)| < 2 


\yg{y)\ dy < eo/2. 


Estimates (|^) and (|6.3|) and Rouche’s theorem imply that A„ G Dn. 

Since the functions <I> and $ assume real values for real A, for every n G M there is 
An between A^ and A^, for which 


(6.5) $(An) = 4>(An) - 4>(An) = (An - A„)<I>'(An). 

We observe that An G Dn and that in view of (|6.2|) and (|6^), 

|$'(An)|>|4>(An)|-£o/2>£o/2. 


On the other hand, 

$(An) 

SO that we get 


$(An) - $(An) = Sn(s') := / S'(l/) siu An?/d|/, 

Jo 

\k-k\D-DP- 
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for all n G N. Recall (see Appendix that the system of functions {sin A„x} forms a 
Riesz basis of H, so that 


n£N 


C > 0 being the corresponding Riesz constant. Combining the above estimates, we 
conclude that 

V ,o . 4C„ . iCe'^ 


nGN 


AnP < ^ 




and the desired continuity follows. 


□ 


Lemma 6.5. The sequence (ak) of norming constants depends continuously on a & 
S+/M. 

Proof. Suppose that a G is in the ^-neighbourhood of a and 5*, are the corre¬ 
sponding norming constants for the operator Tj. As in the proof of Lemma we see 
that 

ll^fc ^/cllll^/c 

where Uk = (/ + Ka)vk and Uk = (4 + K^)vk, and are transformation operators, 
and Vk = v^sinAfco;, Vk = \/2smXkX. Since the norms \\uk\\ are bounded uniformly in 
a G Og, it remains to show that the sequence il-Ufc—Mfc|| is in £2 and > 0 

as e —>• 0. 

We have 

Uk - Uk = Vk - Vk + {K„ - K^)vk + K^{vk - Vk). 

The system {vk) with \k = nk + p.k and (/ifc) G £2 depends continuously on (/Xfc) in the 
sense that 

\\vk -Vkf <2Y \\Pk - Pkf- 

fceN fceN 

Since the norms ||A" 3 ?|| are uniformly bounded in a G O^, also 

Y \\K^{vk - hfe)|P < Cl ^ ll/ifc - ]lkf. 

fceN fceN 

Finally, we observe that the system (vk) is a Riesz basis of TC, so that 


Y\\{Ka-K^)Vkf <C2\\K^-K, 


2 

O' 11 ©2 ’ 


ken 


see Appendix Since K^r depends continuously in ©2 on a G by the results of 
the required continuity follows. □ 


Combining the above lemmata, we arrive at the following result. 


Theorem 6.6. The mapping 

S+/R9d^{(A2),(a,)}G5P 

is a homeomorphism. 

Fixing the spectrum (A^), we can say even more about the corresponding isospectral 
set (cf. [^, Ch. 4] for the regular case g G TY). 

Theorem 6.7. Suppose that the sequence (A^) satisfies assumption (Al). Then the 
set of all isospectral potentials in hF2’"^(0,1) with the Dirichlet spectrum (A^) is analyt¬ 
ically diffeomorphic to an open subset of the Hilbert space £ 2 . The diffeomorphism is 
performed through the sequence {(dk), where j3k = ak — f- 
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Proof. It suffices to note that the correspondence {(5k) 0, 

0(5) = (cOSTt/cS — COsAfcs) — E ( 5 k cos XkS, 
fceN feeN 

where (5k = (5k/ (l+/3fc), is bounded, affine (thus analytic) in {(5k) G £ 2 - Since ak = l+Pk 
are uniformly bounded away from zero, the mapping {(5k) 1 —>■ 0 is analytic in {(5k) G £ 2 , 
and the result follows from Lemmata |6.2| and |6-3i □ 

Analogous stability and isospectrality results hold true for other boundary conditions 
considered in the next section. 

7. The case of other boundary conditions 

The solution of the inverse spectral problem presented in Sections |]-|^ can easily be 
adapted to the case of other types of boundary conditions, e. g., the boundary condi¬ 
tions of the third type at both endpoints, Dirichlet-Neumann, or Neumann-Dirichlet 
ones. Below, we shall briefly discuss the modifications to be made for these cases. 

For cr G 7-f and H,h E C, we consider a Sturm-Liouville operator given by the 

differential expression of ( |1.2|) and the boundary conditions 

mW(0) - Hu{£)) = 0, mW(1) + hu{l) = 0. 

(We recall that = u'{x) — a{x)u{x) is the quasi-derivative of u.) More precisely, 

Ta,H,h is given by 

Ta,H,hU = la{u) := - CTU 

on the domain 

2?( W) = {u G W/{0, 1) I uW G W/{0, 1), IM G n, 

uW(0) - Hu{0) = mW(1) + hu{l) = 0}. 

Observe that Ta^H,h = so that without loss of generality we may (and will) 

assume that H = 0. 

It is known that for all real-valued a E PC and h E M. the operator To-0,/1 is 
selfadjoint, bounded below, and has discrete simple spectrum (A^), k E N. As earlier, 
upon adding a suitable constant to the potential q = a', we can make all the eigenvalues 
positive. Denote by Uk the eigenfunction of To-,o,/i corresponding to the eigenvalue A^ 
and normalized in such a way that Mfc(O) = \/2, and put ak ■= ||Mfc|P- 

Our aim is to solve the inverse spectral problem for To-^ 0 ,/ 1 , i-e., firstly, to describe 
the set SV of all spectral data {(A^), ( 0 ^)} that can be obtained by varying real-valued 
a E Pi and h eM. and, secondly, for given spectral data {(A^), (a^)} G SV, to find the 
corresponding operator To- 0,/1 (i.e., to hnd the corresponding primitive a E Pi oi the 
potential q and the number h G M). 

Lemma 7.1. Suppose that a E Pi is real-valued, h G M, and that Xf < < ... 

are eigenvalues of the operator with A^ > 0. Then Xk = 7r{k — 1) + qik, where 

(/ifc) G £ 2 - 

Proof. Denote by T„q the extension of the operator obtained by omitting the 

boundary condition at the point x = 1. By 0, the operators and are similar, 

and the similarity is performed by the transformation operator I -|- here iLo-,o is 

an integral operator of Volterra type, {Kaflu){x) = kafl{x,y)u{y) dy, and the kernel 
/Co -,0 has the property that k„fl{x, ■) is an 7-f-function for every x E [0, 1]. 
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Put m( •, A) = (/ + Ka^fi)v{ •, A), where v{x, A) = \/2cos Ax. Then Uk = u{-, A^) and 
the numbers ±Xk are solutions of the equation m[^](1,A) + hu{l, X) = 0. Using the 
properties of the transformation operator iP^- o (see Remark 4.3]), we hnd that 

A) = —a/ 2A sin A — \/2A / 5^1 (x) sin Ax dx + / 5f2(a^)\/2 cos Ax dx + C 

Jo Jo 

for some functions gi,g 2 from 7d and a real constant C. It follows that Xk are zeros of 
the analytic function 

<hi(A) := —A sin A + h cos A — A / 5 fi(x) sin Axdx + / 5 f 3 (x) cos Ax dx + U/v^, 

Jo Jo 

where (73 := g 2 +hk(j^Q{l, ■) G 7d. Now the standard analysis (cf. the proof of Lemma |2.3|) 
yields the asymptotics required. □ 

We next establish the asymptotics of the norming constants 

Lemma 7.2. Suppose that a Ehi is real-valued, h G M, and that Uk are eigenfunctions 
of the operator normalized as above. Then Uk = 1 + Pk, where the sequence (Pk) 
belongs to £ 2 - 

Proof of this statement is completely analogous to the proof of Lemma p.4| . The 
minor changes to be made concern notations, namely, we should put Ufc(s) = u(s, A^) = 
A/2cosAfcS and = \/2cos[7r(fc — l)s]. 

Suppose that the sequences (A^) and [ak) are the spectral data for an operator Ta-,o,h 
with a E H and h G M and put 

N ^ 

F := s-hmY^(- {■,Vk)vk - {■,Vk,o)vk,o). 

N^oc-^^Kak / 

k=l 

It is easily seen that / + F is uniformly positive. As in Lemma p.2| , it can be proved 
that F has the following properties. 


Lemma 7.3. The operator F is a Hilbert-Schmidt integral operator with kernel 

(7.1) f{x,y) = (p{x + y) + (p{x-y) 
where 

00 

(7.2) 4 ,( 3 ) = Y. 

k=l 

is an L 2 { 0 , 2 )-function. 


(— cos AfcS — cos 7r(/c — l)s 
\ak 


Next we show that F is related to the transformation operator I + K„q through 
the GLM equation (^Tf) , which can be used to determine uniquely through the 
formula = — (J + Vp)~^V^F G 6^, see Section |[ 

Suppose now that the sequences (A^) and (cr^) consist of positive numbers and are 
as claimed in Lemmata |7.1| and |7.2| . We solve the GLM equation (|4.1|) for K, and then 
in the same manner as in Theorem we show that there exists a E Ti such Uiat 
I + K is the transformation operator for the Sturm-Liouville operators To-,o and Tq^q. 
The only thing that remains to be proved is that our initial sequences (A^) and («*,) 
are the spectral data for the Sturm-Liouville operator Tup^h with some h G M. 

We use identity 


I) and replicate the arguments of the proof of Theorem |5.3| to 
show that none of Uj{l) vanishes as otherwise all Uj{l) would vanish and Xk would have 
a different asymptotics (that of the Neumann-Dirichlet boundary conditions case, see 
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below). Therefore Uj{l) 7 ^ 0 for all j G M, and relation ( f).4| ) implies that there exists 
a constant h G M such that 


(7.3) uf{l)/uj{l) = h 

for all j G N. Thus (A^) are eigenvalues of the Sturm-Liouville operator — ^ + cr' 
subject to the boundary conditions |/f^^( 0 ) = 0 and = hy{l), and the proof is 

complete. 

We summarize the above considerations in the following theorem. 


Theorem 7.4. For two sequences (A^) and {ak) to he the spectral data of a (positive) 
Sturm-Liouville operator Ta-fi^h with real-valued potential q = a' from 14^2“^(0,1) and 
h eM., it is necessary and sufficient that ak are as in (A2) and Xk satisfy the following 
assumption: 


(AT) Xk are all positive, strictly increase with k, and obey the asymptotic relation 
Xk = vr(/c — 1) + /ifc with some i 2 -sequence {fik)T=i- 
The corresponding operator is uniquely recovered from the spectral data through 
formula ( b-3|) , in which the functions 0 and f are given by ( [f.2|) and o respectively, 
and k is the kernel of the solution K of the GLM equation (|T). The number h in the 
boundary conditions is given by ( 0 ). 


In a similar manner the case where one of the boundary conditions is a Dirichlet one 
and the other one is of the third type can be treated. It suffices to consider only the 
cases where H = 00 , h = 0, or H = 0, h = 00 a.s other situations reduce to one of 
these in view of the relation T^^u^h = T„+h',H-h',h+h'- 

The operators To-, 00,0 and To-, 0,00 can be uniquely recovered from the spectral data, the 
sequences of eigenvalues and norming constants. The eigenvalues A^ of the operators 
Ta,oo,o and To-, 0,00 obey the asymptotics Xk = 7r(/c — 1/2) + fik for ^ 2 -sequences (/Xfc), 
while the norming constants ak (defined as in Section |] or Section according as 
H = (yo or H = 0) satisfy (A2). Thus the sets of spectral data for the families of 
Sturm-Liouville operators To-, 00,0 and T^p^ooi where a runs through 74, admit explicit 
descriptions, while the reconstruction algorithm remains the same. The details can 
be recovered by analogy with the analysis of Sections (cf. [|^, Sect. II. 10] for the 
regular case). 


Appendix A. Riesz bases 

In this appendix we gather some well known facts about Riesz bases of sines and 
cosines (see, e.g., [T^ and references therein for a detailed exposition of this topic). 

Recall that a sequence {en)’^=i in a Hilbert space 74 is a Riesz basis if and only if 
any element e E Ti has a unique expansion e = with (c„) G f'2- If (sn) 

is a Riesz basis, then in the above expansion the Fourier coefficients Cn are given by 
Cn = (e, e(j), where is a system biorthogonal to (e^), i.e., which satisfies the 

equalities (efc,e/) = 5kn for all k,n e'H. Moreover, the biorthogonal system (e/) is a 
Riesz basis of 74 as long as (e„) is, in which case for any e G 74 also the expansion 
e = takes place. In particular, if (e„) is a Riesz basis, then for any e G 74 

the sequence (c(^) with c'^ := (e, e„) belongs to £ 2 - 

Also, any Riesz basis (e„) is equivalent to an orthonormal one, i.e., there exists a 
homeomorphism U such that {Uen) is an orthonormal basis of 74. As a result, there 
exists a constant C > 0 (the Riesz constant) such that, for any sequence (ck) E £ 2 , 





in particular, the series Yhck^k converges in Ti for any f' 2 -sequence (cfc). 


Proposition A.l ([]^). Suppose that the real numbers fik, A; G N, tend to zero and 
that the sequence {nk + p,k) strictly increases. Then each of the following systems forms 
a Riesz basis of L2{0, 1); 

(a) {sm{Trkx + 

(b) {sin(7r[A; - l/2]a: + 

(c) {cos(7rA:a; + /ifca;)}^o; 

(d) {cos(7r[/c +l/2]a; + /ifca;)}^o. 


Corollary A.2. Under the assumptions of Proposition pn| the systems (a)-(d) con¬ 
stitute Riesz bases of their closed linear spans in L2(0,2). 


Appendix B. The ideal of Hilbert-Schmidt operators 

In this appendix we recall some necessary facts about the Schatten-von Neumann 
&p ideals (see details in | 0 ). 

Suppose that T is a compact operator in a Hilbert space TC; then |T| = (T*T)^/^ is 
a nonnegative selfadjoint compact operator. Denote by Ai(|T|) > A 2 (|T|) > ... the 
eigenvalues of |T| in non-decreasing order and repeated according to their multiplicity. 
Then Afc(|T|) is called the k-th s-number of T and is denoted by Sk(T). 

The ideal &p, p G [1, 00 ), consists of all compact operators for which the expression 

l|r|le,:= 

is finite. Being endowed with the norm || ■ He^, the ideal Gp becomes a Banach space. 

In particular, Gi is the ideal of trace class operators. For any T G ©1 its matrix trace 
trT := with respect to any orthonormal basis (cfc) coincides with the 

spectral trace, i.e., the sum of all eigenvalues of T repeated according to their algebraic 
multiplicity, and is finite. 

For p = 2 the ideal 62 consists of all Hilbert-Schmidt operators. For T G ©2 and 
any two orthonormal bases {cj) and (e(,) of TL, we have 

ITfa = E l(re,-.4)t = E < oo- 

yfcsN jgn 

If [cj) is a Riesz basis of TL, then {Uej) is an orthonormal basis for some homeomor- 
phism U, and for any T G ©2 the estimate 

ElTejf = l|n7-‘llL<lir|liJT-'lt 

jeN 

holds. Also, ST is a trace class operator whenever S and T are Hilbert-Schmidt ones. 
Moreover, ©2 is a Hilbert space with the scalar product given by 

( 5 ,T )2 := trFT* = trT*R 

If Ti = 1 ^ 2 ( 0 ,1), any Hilbert-Schmidt operator T is an integral one with kernel t given 

by _ 

t{x,y)-= ^{Tej,ek)ek{x)ej{y)] 

here (e^) is any orthonormal basis of TC and the series converges in T 2 (( 0 ,1) x (0,1)). 
Moreover, 

['\tix,y)\^dxdy=\\T\\l^. 

Jo Jo 
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